Draft version February 21, 2012 

Preprint typeset using I^'T^]X style emulatoapj v. 5/2/11 



FRAME-DRAGGING, DISK WARPING, JET PRECESSING, 
AND DIPPED X-RAY LIGHTCURVE OF SW J1644+57 

Wei-Hua Lei^'^, Bing Zhang\ He Gao^ 

Draft version February 21, 2012 

ABSTRACT 

The X-ray transient source Sw J1644-t-57 recently discovered by Swift is believed to be triggered 
by tidal disruption of a star by a rapidly spinning supermassive black hole (SMBH). For such events, 
the outer disk is very likely misaligned with respect to the equatorial plane of the spinning SMBH, 
since the incoming star before disruption most likely has an inclined orbital plane. The tilted disk 
is subject to the Lense-Thirring torque, which tends to twist and warp the disk due to the Bardeen- 
Petterson effect. The inner disk tends to align with the SMBH spin, while the outer region tends 
to remain in the stellar orbital plane, with a transition zone around the Bardeen-Petterson radius. 
The relativistic jet launched via the Blandford-Znajek mechanism from the spinning SMBH would 
undergo precession. The X-ray lightcurve of Sw J1644-I-57 shows a quasi-periodic (2.7-day) variation 
with noticeable narrow dips. We numerically solve a warping disk solution and propose a jet-processing 
model by invoking a Blandford-Znajek jet collimated by a wind launched near the Bardeen-Petterson 
radius. Through simulations, we show that the narrow dips in the X-ray lightcurve can be reproduced 
for a range of geometric configurations. From data we infer that the inclination angle of the initial 
stellar orbit is in the range of 10° — 20° from the SMBH equatorial plane, that the jet should have 
a moderately high Lorentz factor, and that the inclination angle, jet opening angle, and observer's 
viewing angle are such that the duty cycle of the line-of-sight sweeping the jet cone is somewhat less 
than 0.5. 

Subject headings: accretion, accretion disks-black hole physics-magnetic fields 



1. INTRODUCTION 

The discover of the hard X-ray transient event Swift 
J16449.3+573451 ("Sw J1644-H57" hereafter. Burrows et 
al. 2011) by the Swift satellite (Gehrels et al. 2004) 
has stimulated a great interest in studying jets launched 
from tidal disruption events (TDEs). The long variabil- 
ity time scale 5t ~ 100 s (Burrows et al. 2011) and 
its location near the center of a 2: = 0.354 host galaxy 
(Levan et al. 2011) link Sw J1644-I-57 to a super-massive 
black hole (SMBH) (Burrows et al. 2011; Bloom et al. 
2011). The sharp onset and gradual fade-away of X- 
ray flux refer to tidal disruption of a star by a dormant 
SMBH (Bloom et al. 2011; Burrows et al. 2011). The 
super-Eddington X-ray luminosity (Burrows et al. 2011), 
bright radio afterglow (Zauderer et al. 2011), as well as 
a stringent historical X-ray flux upper limit suggest that 
a relativistic jet is launched from a SMBH during the 
TDE, which is not expected in most previous TDE stud- 
ies (e.g. Kobayashi et al. 2004; Strubble & Quataert 
2010; Lodato & Rossi 2011, but see Lu et al. 2008; Gao 
et al. 2010; Giannios & Metzger 2011). ModeUing the 
emission of Sw J1644-I-57 suggests that the jet is highly 
"particle starved" (Burrows et al. 2011), favoring a mag- 
netically launched jet, likely launched via the Blandford- 
Znajek (1977, hereafter BZ) mechanism (Lei & Zhang 
2011), which extracts the spin energy of the BH through 
a magnetic field connecting the BH event horizon and 
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a remote astrophysical load. From observational con- 
straints. Lei & Zhang (2011) found that the SMBH in 
Sw J1644-I-57 carries a moderate to high spin. 

For such events, since the star initially has no knowl- 
edge about the BH spin orientation before being dis- 
rupted and accreted, it is most likely that the initial 
stellar orbit is mis-aligned with the equatorial plane of 
the spinning BH. A natural expectation is that at least 
the outer part of the accretion disk is also misaligned. 
The tilted disk surrounding a spinning BH is subject 
to the Lense-Thirring (Lense & Thirring 1918, LT here- 
after) torque. The combined action of the LT torque and 
the internal viscosity of the accretion disk would lead to 
a twisted and warped disk, with the inner part of the 
disk bent towards the BH equatorial plane due to the 
frame-dragging effect. This is known as the Bardeen- 
Petterson effect (Bardeen & Petterson 1975, hereafter 
BP). For a fully developed BP disk, the inner part of the 
disk tends to be aligned with the BH equatorial plane, 
while the outer part of the disk remains aligned with the 
original stellar orbital plane. The transition radius be- 
tween these two regimes is the BP radius. The warped 
disks have been directly observed by water maser obser- 
vations in NGC 4258 (Miyoshi et al. 1995; Neufeld & 
Maloney 1995; Herrnstein, Greenhill & Moran 1996) and 
the Circinius galaxy (Greenhill et al. 2003). The ap- 
parent lack of correlation between the direction of the 
radio jets emanating from active galactic nuclei (AGN) 
and the plane of the host galaxy disks (Kinney et al. 
2000; Schmitt et al. 2002) can also be explained by disk 
warping. 

Besides this theoretical motivation, some curious ob- 
servational facts also point towards the possibility of a 
warped disk. Burrows et al. (2011) reported a rough 
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period of 230 ks 2.7 days) in the long-term X-ray 
lightcurve, although with a < 3(t significance. The sim- 
ilar rough period was also revealed recently by Saxton 
et al. (20f 2). More importantly, a closer inspection sug- 
gests that the late X-ray lightcurve show interesting dips. 
A natural mechanism to account for these dips would be 
a precessing jet, possibly caused by a warped BP disk. 
Stone & Loeb (2012) recently also discussed the similar 
topic, but didn't pay attention to these peculiar observa- 
tional effects. They led to the conclusion that the stellar 
orbital plane is almost aligned with the BH equatorial 
plane due to the lack of observational evidence of jet 
processing. 

In this paper, we develop a model to interpret the 
quasi-periodic feature and the narrow dips observed in 
the lightcurve of Sw J1644-I-57 within the framework of 
the BP mechanism. In Section 2, we first apply a math- 
matical method, i.e. Stepwise Filter Correlation (SFC) 
method (Gao et al. 2012), to study the X-ray lightcurve 
data of Sw 1644-1-57. We confirm that a ~ 2.7-day quasi- 
periodic signal exists in the late lightcurve after the ini- 
tial flaring phase. We then develop a numerical model in 
Section 3 to study disk warping and jet processing. In 
Section 4, we apply the observational data to constrain 
the model parameters, and found several parameter sets 
that can roughly reproduce the observations. Our con- 
clusions are presented in Section 5. 

2. OBSERVATIONAL EVIDENCE OF JET PRECESSION 

Even though not with high significance, two pieces of 
evidence hint towards the possibility of a precessing jet 
in Sw J1644-I-57. The first evidence is the noticeable dips 
in the XRT lightcurve at late times. The upper panel of 
FigH] shows the XRT lightcurve of Sw J1644-h57 after 5 
days since the first trigger. Some deep dips are clearly 
visible. The dashed vertical lines mark the 2.7-day period 
(see more below) . One can see that the dips roughly track 
this period. Dips are not easy to explain within models 
invoking emission, but would be naturally interpreted 
within the frame work of jet precession where the line of 
sight is allowed to move away from the emission beam. 

The second piece of evidence is the rough 2.7-day peri- 
odicity, which was pointed out by Burrows et al. (2011) 
and confirmed by Saxton et al. (2012). We investigate 
this independently using a new method. Stepwise Fil- 
ter Correlation (SFC) method, recently developed in our 
group (Gao et al. 2012). The SFC method is a sig- 
nal processing algorithm which stepwise filters signals 
above a frequency, and looks for correlation between the 
lightcurves of adjacent filterings. If there is a spectral 
component around a particular frequency, the lightcurves 
before and after filtering this particular frequency can 
be very different, leading to a dip in the correlation 
curve. The detail of this method and its application 
to GRB lightcurves are presented in Gao et al. (2012). 
Here we apply the method to the XRT lightcurve of Sw 
J1644-I-57. We first perform an analysis to the entire 
lightcurve (curve (a) of lower panel of FiglT|). Since early 
lightcurve displays many erratic, bright X-ray flares, the 
signal is dominated by these high-frequency components, 
and the 2.7-day component cannot be identified. These 
flares may be related to fallback of the stellar debris with 
high and fluctuant accretion rate (Cannizzo et al. 2011). 
If we remove the early data (say, the first two days) , the 



2.7-day dip clearly shows up (curve (b) of lower panel 
of Fig[T]). We then gradually remove more and more 
early data points and everytime redo the SFC analysis 
(curves (c-e) of the lower panel of FigH]). The 2.7-day 
dip is always there. Our analysis suggests that, consis- 
tent with the results of Burrows et al. (2011) and Sax- 
ton et al. (2012), that the late time X-ray lightcurve of 
Sw J1644-I-57 has a quasi-periodic signal with a period 
roughly 2.7 days. 
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Fig. 1.— Upper panel: Swift XRT light curve of Sw J1644+57 
from 5-30 days, which show clear dips. The dashed line shows 
the 2.7-day period. Data is from Burrows et al. (2011). Bottom 
panel: the stepwise filter correlation (SFC) analysis of the X-ray 
lightucrve following Gao et al. (2012). The curve (a) is the result 
for the entire lightcurve, curve (b), (c), (d), (e) and (f) are the 
analyses by manually removing the first 2, 5, 10, 15, and 20 days 
of the data, respectively. A dip corresponding to ~ 2.7 days is 
evident once the first 2-day of data is removed. 

3. WARPED ACCRETION DISK AND JET PRECESSION 

In simple physical terms, disruption of a star occurs 
when it comes to a SMBH closer than the tidal disrup- 
tion radius Rt, which is determined by demanding that 
the mean density of space volume enclosed by Rt, i.e., 
M,/(47ri?^/3), is equal to the density of the star, where 
M, is the mass of the BH. The tidal disruption radius is 
then given by 

i?T i^y^'R. ^ 7 X lO^W'/'nAfiZ/cm 

~ 47i?gTO^^/^r*M.~g/^ (1) 
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Fig. 2. — The sketch of the model. A star is tidally disrupted as it 
gets close to a dormant SMBH. The star orbit is generally inclined 
with the SMBH spin vector Jbh- After disruption, an accretion 
disk form, and is warped due to the Bardeen-Petterson effect. The 
jet is expected to precess with an angle 9; around the BH spin 
axis. For a jet with half opening angle Oj and an observer located 
at Ooha from the BH spin axis, the angle ip between the observer's 



line of sight 



and jet axis fj^t is changing with time. The line 



of sight can exit llic jet cone at certain phases, leading to dips in 
the lightcurve. The Bardeen-Peterson radius R^p is marked. The 
jet direction is defined by the disk normal direction around R^p . 

where = M^/Mq, n = R*/Rq, M,^q = M./IO^M© 
and i?g = GM,/c^. For a solar-likc star (M, = Mq 
and i?, = i?©) disrupted by a 2 x IQ^Mq BH, the tidal 
disruption radius is approximately 30i?g from the BH. 

Stellar disruption occurs if R^s < -Rp < -Rt, where i?p 
is the periastron distance from the star, and iZms is the 
marginally stable orbit (Bardccn ct al. 1972). The star 
should be swallowed whole instead of being disrupted if 
it goes inside i?ms- We define n = Rp/Rt- 

We idealize the disk as extending from iJms to the cir- 
cularization radius -Rout — vRp^ where 77 ~ 2. In our cal- 
culation, the typical values of the parameters are r? = 2, 
= 1, = l,and n = I. 

The dragging of the incrtial frame (frame dragging) 
produced by a Kerr black hole causes precession of a 
particle if its orbital plane is inclined with respect to the 
equatorial plane of the black hole. This effect is known 
as Lense-Thirring precession. The precession angular ve- 
locity fluT is given by (e.g. Wilkins 1972) 



^LT(i?) - 



(2) 



where J, = a,GM^/c\s the BH angular momentum, and 
a, is the BH spin parameter. 



The precession period Tp can be estimated as 



27T/n 



LT 



(3) 



For a BH with mass 2 x IO^Mq and spin a, = 0.9, we 
have Tp ~ 2.7 days at iZ = 19iZg. 

When an accretion disk docs not align with the equato- 
rial plane of the BH, i.e., when the angular momentum of 
the accretion disk L is misaligned with respect to the di- 
rection of J,, the LT effect causes a precession effect that 
twists the disk plane due to the coupling of J, with L. 



The torque tends to align the angular momentum of the 
matter in the disk with that of the BH, thus causing the 
inclination angle between the angular momentum vectors 
to decrease with decreasing distance from the BH. This 
is the BP effect (Bardeen & Petterson 1975), which is the 
combination of the LT effect and the internal viscosity of 
the accretion disk. The effect is more prominant in the 
inner part of the disk due to the short range of the LT 
effect. The outer part of the disk tends to remain in the 
orientation defined by the original stellar orbit. 

Considering that the disk may be thick. Stone & Loeb 
(2012) suggested that the disk is not subject to the BP 
warping effect. This conclusion is based on the GRMHD 
simulations by Fragile et al. (2007), which invoked a 
weak magnetic field to account for angular momentum 
transport. The dynamical viscosity and magnetic stress 
in the disk are subject to large uncertainty, which can 
be larger than invoked by Fragile et al. (2007). For 
relativistic twisted disk around a Kerr BH, Zhuravlcv & 
Ivanov (2008) found that BP effect can take place if the 
viscosity parameter a is sufficiently large. Motivated by 
the observational evidence discussed in Section 2, in this 
paper we suggest that the BP effect indeed plays a role, at 
least after the initial flares when the accretion rate drops 
significantly. To prodticc the high luminosity at early 
phase, the accretion rate is around lO^''Af0.s^^ (Lei & 
Zhang 2011; Shao et al. 2011). After this phase, the 
accretion rate declines to ~ 1O~^M0S~^, about 2 orders 
lower. For a BH with mass IO^Mq and the disk efficiency 
about 10%, the Eddington accretion rate is around ~ 
lO^*M0.s^^. Therefore we could assume that the disk 
becomes thin after the flare phase, and the BP effect 
would come into play to warp the disk. 

For a warped disk, we follow the formalism of Pringle 
(1992). There are two viscosity parameters and 1^2, 
where vi is the standard shear viscosity in a flat disk, 
and f2 is the viscosity associated with the vertical shear 
motion describing the diffusion of warping distortion 
through the disk. In order to study the BP effect, We re- 
write Pringle's equation by adding an effective coupling 
due to the LT precession, fli/r x L (Scheuer & Feiler 
1996; Lodato & Pringle 2006): 



dL 

'dt 



3 d Q 



RdR' 
1 d 



S dR 



{v^^R^'^)^ 



d\ 



1 d 



1 dl 



(4) 



where L = |L| = Y^^GM^R is the angular momentum 
per unit area of disk, S is the surface density, 1 is the unit 
vector indicating the local direction of the speciflc angu- 
lar momentum in the disk. Note that i? is a spherical 
coordinate, and is not the cylindrical radius. 

The transition radius between the two regimes is 
known as the Bardeen-Petterson radius i?BPj which is 
the radius at which the warping propagation time-scale 
= R^/1^2 equals the local forced precession rate Tp, 
giving (Scheuer & Feiler 1996; Caproni et al. 2007) 



Rbp = 



V2 



^\it{R^v) 



(5) 
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Here again V2 is the viscosity normal to the accretion 
disk. Roughly speaking, the disk becomes aligned with 
the BH spin at i? <C -Rbp and keeps its original inclina- 
tion for R Rbp ■ 

Scheuer & Feiler (1996) have found an analytical 
steady-state solution to the above equation (removing 
the time-dependent term dh/dt), subject to various sim- 
plifying assumptions about the disk density distribution 
and the inner boundary condition. This solution may be 
inaccurate at very small radii from the BH. It is based 
on a first-order approximation and a small inclination 
angle (^orbit 1) assumption. In our modeling, instead 
of applying this solution, we solve the equations numeri- 
cally following the numerical scheme outlined in Pringle 
(1992). 

The relation between vi and z/2 is the main uncertainty 
of the problem. For simplicity, we take vxjvi — 1 in the 
calculation. The solution for the inclination angle at time 
t and radius R depends on the BP radius i?BP , the outer 
disk boundary i?out and the orbital inclination 6'orbit- In 
order to set up a twisted disk, we add matter with an in- 
clination angle 6'orbit to the grid near the outer boundary 
^out- We run the code for about 10^ time-steps without 
the LT effect to settle the disk into a numerically equi- 
librium state. We then evolve the disk by applying the 
forced precession due to the LT torque. The disk twist- 
ing evolution slows down rapidly. After about 2 x 10^ 
runs, the inclination no longer evolve significantly, and 
we take this as the steady state. In FigEl the numerical 
solutions in the steady state are shown for different i?BP, 

-Rout and ^orbit- 
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Fig. 3. — Disk inclination fiinc vs. radius _R for different -RbPi 
Rout and ©orbit: Thick solid line: -Rbp = 10-Rg, -Rout = 50-Rg 
and ©orbit = 30°; thin solid line: -Rbp = 10-Rg, -Rout = 60-Rg 
and ©orbit = 30°; dashed line: -Rpp = 25-Rg, -Rout = 60-Rg and 
^orbit = 30°; and dotted line: -Rbp = 25-Rg, -Rout = 60-Rg and 
f^orbit = 20°. 

From FiglSl we find that the inclination of the inner 
disk critically depends on the BP radius -Rbp and the 
initial inclination of stellar orbit 6'orbit- The effect of the 
outer disk boundary is not significant for the inner disk. 

A precessing disk would induce jet precession. In our 
model, the relativistic jet is launched near the BH by 
the BZ process (Lei & Zhang 2011). However, the jet di- 



rection is not determined by the normal direction at the 
inner disk, but is rather defined by normal direction at 
an outer region in the disk near R^-p . This is because the 
magnetic field threading the disk would drive an outfiow 
from the disk surface via the Blandford & Payne (1982) 
mechanism. A centrifugally driven outfiow is possible if 
the poloidal component of the magnetic field makes an 
angle less than 60° with the disk surface. This is most 
likely to occur at -Rgp, where the disk inclination angle 
significantly changes. A similar proposal was introduced 
to study other BH accretion systems (e.g. Begelman et 
al. 2006) and gamma-ray bursts (e.g. Lei et al. 2007). 
In the following calculation, we suggest that the preces- 
sion period and inclination angle of the relativistic jet is 
defined by the physical conditions at -Rbp- 

4. OBSERVATIONAL CONSTRAINTS 

Some constraints on the model can be derived from the 
observational data outlined in Section 2. 

The first constraint is the 2.7-day quasi-period. In the 
precession model delineated above, this period should 
correspond to the LT period at the BP radius, -Rbp. Ac- 
cording to Lei & Zhang (2011), the most probable value 
of the BH spin in Sw J1644-h57 is a, = 0.9. We take 
this value for the BH spin in the calculation. Applying 
Tp = 2.7days into Eq.®, we can estimate the BP radius 
-Rbp that satisfy the observations (FigU]). For compari- 
son, we also show the disk outer edge -Rout in FiglH where 
M, = Mq, -Rh. — Rq, k = 1 and 77 = 2 are adopted as 
the typical values. It is found that -Rbp is located close 
to the BH. 
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Fig. 4. — The constraint on the Bardeen-Petterson radius -Rbp 
(solid line) as a function of BH mass, based on the requrement of 
rp{-RBp) = 2.7days with a, = 0.9. The disk outer edge -Rout is 
also shown as the dashed line. 

The next observational constraint comes from the ob- 
served shape and depth of the dips seen in the lightcurve. 
They are defined by the complicated combinations of sev- 
eral unknown parameters, including the bulk Lorcntz fac- 
tor r, the jet half opening angle 0j, the inclination angle 
of the jet 6'i, as well as the observer's viewing angle 6'obs 
with respect to the BH spin axis. In the following we 
discuss these constraints in detail. 



Jet precessing in Sw J1644+57 



5 



At any time t, the angle between the observer fobs and 
jet axis fjot(0 (the axis of the presumed conical jet) is 
denoted by ^{t) (see FiglU, which is defined as 



(6) 



A dip in the lightcurve can be seen if the line of sight is 
outside the jet cone, i.e., tlj{t) > 9j. 

For simplicity, we start with a uniform conical jet, so 
that when ip < 9j the observer receives a uniform flux. 
For an off-axis observer at ■0 > ^jj the observed flux 
density can be written as (Granot et al. 2002) 



F,{,p,t)^D^F,/D{0,Dt), 



(7) 



where D = {1 - I3)/{1 - I3cos0), (3 = ^1 - l/r2 and 
ij) = max(0,^obs — ^'j)- Define F{ip,t) — vF^{'ij),t), we 
have the X-ray flux at ip 



F{i;,t) = D^F{0,Dt). 



(8) 



For a tidal disruption event, the long-term lightcurve is 
defined by the fall-back accretion rate, which is oc i"^/"^. 
The long term X-ray lightcurve of Sw J1644-I-57 is indeed 
consistent with this behavior, i.e. F{t) cx t~^/^. We 
therefore have 



F{'iP,t) = D^^'^F{Q,t). 



(9) 



The shape the lightcurve depends on several factors. 
The first one is the Lorentz factor. In FigIS] we show 
how the Lorentz factor affects the observed structure of 
jet. Even if the physical jet has a sharp cutoff near the 
edge, the observed jet structure has a wing defined by 
the gradual decrease of the Doppler factor. If the jet 
Lorentz factor is large, the flux drops sharply as the line 
of sight goes away from the jet cone, and one would see a 
very deep dip. On the other hand, the dip would be very 
shallow if the jet has a smaller Lorentz factor. In the 
observed lightcurve, the dips do not have uniform depth, 
and in some periods the dips disappear completely. This 
might be caused by fluctuation of the jet Lorentz factor 
as a function of time. 




Fig. 5. — The flux as a function of view angle for the jet for 
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The second parameter to shape the peaks and dips in 
the lightcurve is the duty cycle 



X = tohs/Tp < I, 



(10) 



where tohs is the time when the observer's line of sight 
is inside the jet cone during one precession period. A 
large value of x would correspond to narrow dips, while 
a small x would correspond to wide dips. This ratio can 
be obtained by solving the equation ip{t) = 6y In general 
one obtains x = x{G\,0ohs,6i)- 
The last parameter 



A(6'i, 6'obs, ^'j, F) 



m 



m 



<i, (11) 



defines how deep the dip is. Besides the three angle pa- 
rameters 6'obs, 9]), the Lorentz factor F also enters the 
problem. 

Finally, the jet opening angle 0j and the mean Lorentz 
factor F of Sw J1644-|-57-like events can be constrained 
by through a statistical argument. Based on the com- 
parison between observed event rate of these events and 
the more general TDE event rate constrained from both 
TDE theory and observational data, one can constrain 
the beaming factor /b (Burrows et al. 2011; Lei & Zhang 
2011) 



/b 



This requires that F > 9.5 and <Q° . In the following 
calculations, we take 9-^ = 6° and F = 20 as the typi- 
cal values. A relative large value of F is demanded to 
reproduce the depth of the observed dips. 

The values of 9i and 0obs can be then obtained by com- 
bining Eqs. (|10|) and (jlip . if we obtain information about 
the values of x and A. In FiglBl the relation between 9i 
and 9ohs for different x and A are shown. Based on ob- 
servations, the value of x in Sw J1644-I-57 cannot be too 
large (e.g., point "d" in Figl6]), since otherwise the dips 
would be too shallow to be observed (e.g. FiglTJi). Sim- 
ilarly, A cannot be too small (e.g., point "a" in Fig (5]), 
since otherwise the peaks and dips would be too sharp 
(e.g. Fig[7^). Observationally, moderate values of both 
parameters, e.g. x 0.4 and A ~ 0.3, are favored (e.g. 
points "b" or "c" in Fig|6l see FigslTjj and [7]:). Since we 
observe most of the jet beam, the line of sight cannot be 
too far away from the BH spin axis, i.e. ^obs is expected 
to be small. From FiglHl these constraints suggest that 
the jet inclination angle 9i is likely similar to ~ 0j. 

Based on the above discussion, the BP radius i?BP 
could be obtained with the observed jet precession pe- 
riod (2.7 days). Putting all the constrains into the nu- 
merical model described in Section 2, we can infer the 
inclination of the stellar orbit before disruption ^orbit- 
The results are shown in Fig[51 From Figl5]we find that 
the inclination angle of the initial stellar orbit can be as 
large as 20°. For small BH masses, ^orbit can be much 
greater than 9i. However, for large BH masses, ^orbit is 
not much larger than 9i, mainly due to a small value of 
R^/Rg involved. 

Finally, the observed X-ray lightcurve Sw J1644-I-57 
does not have strict periodic dips. In some phases the 
dips disappear, and in some phases the dips do not ap- 
pear at exactly the predicted phase (Burrows et al. 2011, 



Robs 
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Fig. 6. — The relation between Bi and 9chs for different x (thin 
lines, solid: x = 0-3; dashed: x = OA; dotted: x = 0-5; dot-dashed: 
X = 0.7) and A (thick lines, solid: A = 0.5; dashed: A = 0.3; dotted: 
A = 0.1), where the jet opening angle 8j = 6° and the Lorentz 
factor r = 20. 

and FigH]). We believe that this is related to some 
stochastic processes involved in accretion and jet launch- 
ing. The Lorentz factor T may vary with time, and R^p 
may slightly vary at different epochs. In order to test 
this idea, we carry out a range of simulations to the 
lightcurves. First, we introduce a small time-scale vari- 
ability in the lightcurve overlapped on the F(0, t) oc t~^^^ 
envelope (Burrows et al. 2011), which may be a result 
of stochastic magnetic dissipation (e.g. Zhang & Yan 
2011) or fluctuations of accretion rate (e.g. Wang et al. 
2006). Next, we introduce a variation of Lorentz factor 
by allowing F to randomly vary in the range of (2, 25) 
during each precession period. We simulate four cases 
in Figl?! which correspond to "a", "b", "c" and "d" in 
FiglHl We find that the change of Lorentz factor can in- 
deed result in the disappearance of dips at some epochs. 
The peaks and dips in FiglTK are well separated and are 
too spiky. The lightcurve in FiglTJl on the other hand is 
too smooth, without clear evidence of dips. This suggests 
that they do not resemble the geometric configuration of 
Sw J1644-f57. The other two cases, i.e., FigsEb and[7h, 
give interesting quasi-periodic features with noticeable 
dips. Even though we cannot fully constrain the geome- 
try of the Sw J1644-f57 system, we can say that it should 
resemble the conditions of these two cases based on the 
simulation results. 

5. CONCLUSIONS AND DISCUSSION 

In this paper, we propose that the Sw J1644-I-57 jet 
launched from a SMBH precesses. This comes from two 
independent arguments. First, a tidally disrupted star 
should not have knowledge about the spin direction of 
the BH before disruption, so that the original stellar orbit 
very likely has an inclination with respect to the equa- 
torial plane of the BH. A rapidly spinning BH (which is 
probably the main reason to launch a relativistic jet. Lei 
& Zhang 2011) tends to distort the accretion flow and 
warp the accretion disk due to strong LT frame drag- 
ging. This would naturally result in a precessing jet if 
the accretion rate is low enough. Second, observationally 
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Fig. 7.— Several simulated XRT lightcurves of Sw J1644+57. 
The four cases correspond to the "a,b,c,d" points of the parameter 
space in Fig. |6] (a): 0obs = 5.0 deg, 0] = 7.4 deg; (b): ^obs = 4.0 
deg, 01 = 6.0 deg; (c): S^bB = 2.1 deg, 5; = 6.3 deg; (d): Sobs = 3.3 
deg, 9i = 3.3 deg. It is shown that the case (a) has too sharp dips, 
while the case (d) essentially smeared out the dips. The cases (b) 
and (c) would be close to what is observed from the source. 



Jet precessing in Sw J1644+57 



7 



20 




10 



2 4 6 8 10 

Fig. 8. — The inferred inclination of stellar orbit (?orbit vs. the 
BH mass, where Tp = 2.7 days. The solid, dashed and dotted lines 
correspond to Si = 7°, 6° and 5°, respectively. 



two pieces of evidence hint towards a precessing jet: the 
rough 2.7-day periodicity and the noticeable lightcurve 
dips. We performed a SFC analysis to the observed X- 
ray lightcurve, and confirmed the 2.7-day rough period. 
We then carried out detailed theoretical modeling on disk 
warping and jet precessing, and inferred the underlying 
parameters of the jet system from the observational data. 
We found that due to disk warping, the original stellar 
orbit does not have to be nearly aligned with the BH 
equitorial plane (cf. Stone & Loeb 2012). Rather, 6'orbit 
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can be as large as 20°. In order to repduce the depths 
of the lightcurve dips and their distribution, the Lorentz 
factor should range from moderate values to values as 
high as 20. The duty cycle x, which describes the frac- 
tion of time when the line of sight intersects the jet, is 
found to be a moderate value, slightly smaller than 0.5. 

Through simulations, we were able to reproduce 
lightcurves similar to the observed one. Nixon & King 
(2012) interprets the diminishing dips in some periods 
as due to jet nutation. Our SFC analysis does not show 
another quasi frequency that might be related to nuta- 
tion. We attribute the diminishing dips to fluctuations 
in r. The mis-match between the 2.7-day perdiod and 
some dips is also understandable, since the outflow that 
collimates the jet may not always launched strictly from 
Rbp, or i?BP could vary stochastically during the event. 

The main uncertainty of our calculation is the viscosity 
parameters i>i and 1^2 ■ We have assumed a constant value 
for both throughout the event. More complicated evolu- 
tion is possible (e.g. Chen et al. 2009). Some efforts to 
constrain these parameters have been made (e.g. Nixon 
& King 2012), and future detailed numerical simulations 
are needed to better understand the viscosity. 
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